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Abstract 

One considers the superfluid (SF) state of a Bose liquid with a strong repulsion between 
bosons, in which at T = 0, along with a weak single-particle Bose-Einstein condensate (BEG), 
there exists an intensive pair coherent condensate (PCC), analogous to the Cooper conden- 
sate in a Fermi liquid with an attraction between the fermions. Such a PCC emerges in a 
system of bosons due to an oscillating sign-changing momentum dependence of the Fourier 
component of the pair interaction potential, which is characteristic of a certain family of 
repulsion potentials, for example, for a regularized "hard spheres" model or for potentials 
with finite jumps or inflection points with an inflnite derivative. In such cases, the Fourier 
component is negative in some domain of nonzero momentum transfer, which corresponds to 
an effective attraction of a quantum mechanical (diffraction) nature. The collective effects of 
renormalization ("screening") of the initial interaction, which are described by the bosonic 
polarization operator, due to its negative sign on the "mass shell", lead to a suppression of 
the repulsion and an enhancement of the effective attraction in the respective domains of 
momentum space. In the process of building a self-consistent model of the SF state, it is the 
ratio of the BEC density to the full density of the liquid no/n <C 1 that is used as a small 
parameter — unlike in the Bogolyubov theory for a quasi-ideal Bose gas, in which the small 
parameter is the ratio of the number of supracondensate excitations to the number of particles 
in an intensive BEC, (n — no)/no ^ 1. A closed system of nonlinear integral equations for the 
normal Eii(p,cij) and anomalous Ei2(p,a;) self-energy parts is obtained, in the framework of 
a renormalized perturbation theory built on combined hydrodynamic (at p — > 0) and fleld (at 
p 7^ 0) variables, whose usage eliminates infrared divergencies and ensures that the functions 
Eij(p, e) are analytic at p and e and that the SF order parameter Ei2(0, 0) 7^ at 
r = 0. 

In the framework of the hard-spheres model, a spectrum of quasiparticles is obtained, 
which is in good accordance with the experimental spectrum of elementary excitations in 
superfluid ''He. It is shown that the roton minimum in the quasiparticle spectrum is directly 
associated with the first negative minimum of the Fourier component of the renormalized 
potential of pair interaction between the bosons. Finally, the question of applicability of the 
Landau criterion to the description of the SF state of *He in the absence of quantum vortices 
is discussed. 

PACS: 67.57.-z 



*On leave of absence. 
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1 Introduction 



Despite the big progress achieved in the theory of superfluidity since the pioneer 
works by Landau Bogolyubov [Q, Feynman and others [Q-||TD[, the task of 



constructing a microscopic theory of a superfluid (SF) state of a ''He Bose hquid 
cannot be considered complete. In fact, such questions of principle as: (i) the origin 
of the roton minimum in the spectrum of elementary excitations; (ii) applicability 
of the Landau criterion of superfluidity for the determination of the critical velocity 
of dissipationless flow in SF helium (He II); (iii) the quantum mechanical structure 
of the SF component of the ''He Bose liquid below the A point, at T < Ta = 2.17 K, 
etc., remain unsolved. In particular, according to the latest results in quantum 
evaporation of ^He atoms the maximal density po of the single-particle Bose- 
Einstein condensate (BEC) in the ^He Bose liquid even at very low temperatures 
T <ti Tx does not exceed 10% of the total density p of liquid ^He, whereas the 
density of the SF component — >• p at T — 0. Such a low density of the BEC is 
implied by a strong interaction between ^He atoms and is an indication of the fact 
that such an "exhausted" BEC cannot by itself form the microscopic basis of the 
SF component p^. Therefore, the quantum structure of the SF condensate in He II 
with the "excess" density (p^ — po) ^ po calls for a more thorough investigation 

0-0- 

On the other hand, numerous precise experiments on the restoration of the dy- 
namic structure factor ^(p, e) in liquid ^He, involving inelastic neutron scattering 
T^-|T^, show that the temperature dependence of the spectrum of elementary 



excitations E{p), associated with collective density oscillations in the ^He Bose liq- 
uid, is very weak all the way up to the A point, at all momenta, including the 
phonon, maxon, and roton bands. This means that the critical velocity, determined 
according to the Landau superfluidity criterion, Vc = min [E{p)/p], hardly changes 
as T goes up and does not tend to zero as T Ta. At the same time, the break- 



down of superfluidity in macroscopic He II flows is known ||T9[ to be implied by 
the processes of creation of Onsager- Feynman quantum vortices or Anderson closed 
vortex rings. As a result, the threshold velocity v* of breakdown of dissipationless 
flow as observed in He II, may be two orders of magnitude less than the critical 
velocity Vc — [Ar/pr] — 60 m/s associated with the roton gap — 8.6 K in the 

quasiparticle spectrum E{p) at the point p = pr — 1-9 A . 

However, under the conditions where creation and motion of vortices (or vortex 
rings) is hindered, much higher values of threshold velocity can be achieved. For 
example, in ultrathin films and capillaries at T < 1 K, maximal values of v* ~ 
(2—3) m/ s were observed [|l^], and for the passage of He II through narrow apertures 
in thin partitions, critical velocities ~ (8 — 10) m/s were registered pO[, pi 



Moreover, in experiments on acceleration of ions in He II at pressures P ~ (15 — 
20) bar, threshold velocities of more than 50 m/s were achieved, close to the roton 
limit Vc — Aj./pr- In this connection, the problem of ab initio theoretical description 
of the roton minimum in the quasiparticle spectrum E{p) in liquid helium, remains 
topical. 



For the first time, such a problem was considered in Ref. p3[ (see also p3]) within 



the "hard-spheres" model, in which the Fourier component of the regularized pair 
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interaction for the S scattering is a sign-changing function of momentum trans- 
fer. Substituting such a potential into the Bogolyubov quasiparticle spectrum of a 
weakly nonideal Bose gas leads, under certain conditions, to the appearance of a 
minimum, analogous to the roton minimum in the empirical spectrum of liquid ^He. 
A similar problem was considered in Ref. [EH] for a "semitransparent spheres" model 
with a finite potential jump. However, for a Bose liquid with strong interaction be- 
tween particles and a suppressed BEC {uq <^ n) the Bogolyubov approximation is 
not applicable. 

The questions discussed in this paper are both those of the quantum structure 
of the SF state of a Bose liquid and the calculation of the spectrum of elementary 
excitations based on a specific form of pair interaction between the bosons. 

Our approach is based on the microscopic model ^6] of superfluidity of a Bose 
liquid with a suppressed BEC and an intensive pair coherent condensate (PCC), 
which can arise from a sufficiently strong effective attraction between bosons in some 
domains of momentum space (see below) and is analogous to the Cooper condensate 
in a Fermi liquid with attraction between fermions near the Fermi surface As 



a small parameter, one uses the ratio of the BEC density to the total Bose liquid 
density {no/n) <^ 1, unlike in the Bogolyubov theory for a quasi-ideal Bose 
gas, in which the small parameter is the ratio of the number of supracondensate 
excitations to the number of density in the intensive BEC, (n — no) /no <^ 1. 

Because of this, the SF state within the model at hand can be described by 
a "truncated" self-consistent system of Dyson-Belyaev equations for the normal 
and anomalous single-particle Green functions Gij{k,uj) and the self-energy parts 
Tiij{k, uo) without account for the diagrams of second and higher orders in the BEC 
density. The renormalized field perturbation theory is used; it is built on 



combined field variables [Q, [^, which in the long- wave limit (p 0) reduce 
to the hydrodynamic variables of macroscopic quantum (at T = 0) or two-liquid 
(at T 7^ 0) hydrodynamics, whereas in the short-wave band they correspond to 
the bosonic quasiparticle creation and annihilation operators. In this case, the SF 
component ps is a superposition of the "exhausted" single-particle BEC and an 
intensive "Cooperlike" PCC with coinciding phases (signs) of the corresponding 
order parameters. 

The pair interaction between bosons was chosen in the form of a regularized 



repulsion potential in the "hard spheres" model p3|, p4|, whose Fourier component 



V{p) is an oscillating sign-changing function of momentum transfer p due to mutual 
quantum diffraction of particles. The same oscillating sign-changing behavior is 
characteristic of Fourier components of potentials with finite jumps or infiection 
points with an infinite derivative. 

As a result of renormalization ("screening") of the initial interaction V{p) due to 
multiparticle collective correlations, which are described by the boson polarization 
operator n(p, cj), the interaction gets suppressed in the domains of momentum space 
where V{p) > and enhanced where V{p) < 0. Such a suppression of repulsion 
and enhancement of attraction is implied by the negative sign of the real part of 
n(p, cj) on the "mass shell" u = E{p) for a decayless quasiparticle spectrum. It is 
shown that the integral contribution of the domains of effective attraction in the 
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renormalized sign-changing interaction 

V{p) = V{p) [1 - V{p)Re (n(p, E{p))]-' 

can be sufficient for the formation of an intensive bosonic PCC in momentum space 
(ahhough not for the formation of bound boson pairs in real space). 

Self-consistent numerical calculations of the boson self-energy and polarization 
operator, pair order parameter, and quasiparticle spectrum at T = 0, involving an 
iteration scheme, have allowed us to find conditions for the theoretical spectrum 
E{p) to coincide with the experimentally observed elementary excitation spectrum 
in ^He. At the same time it is shown that the roton minimum in the quasiparticle 
spectrum E{p) of a Bose liquid is directly associated with the first negative minimum 
of the Fourier component of the renormalized potential of pair interaction between 
the bosons (akin to the minimum in the Bogolyubov spectrum of a weakly 
nonideal dilute Bose gas |]23[|-||26[|). 



Finally, the question of applicability of the Landau criterion to the description 
of the SF state of ^He in the absence of quantum vortices is discussed. 



2 Green functions and equations for the self-energy parts 
in the model of a Bose liquid with a suppressed BEC in 
the renormalized perturbation theory 

The main difficulty of the microscopic description of the SF state of a Bose liquid 
with a nonzero BEC is the fact that applying the renormalized perturbation theory 
directly |^ leads, as was shown in [|T^- |T^, to a whole number of divergences at 
small energies e — > and momenta k ^ and, as a consequence, to erroneous 
results in the calculations of different physical quantities. 

Thus, for example, for a Bose system with weak interaction, when the ratio of 
the mean potential energy ^(^0)^0 (^0 being a typical momentum transfer) to the 
corresponding kinetic energy k'^/2m of the bosons is small, the zeroth-approximation 
polarization operator n(k, cu) and the density-density response function Il(k, cj) 
calculated to the ffist order in the small parameter of interaction ^ = mkoV{kQ) <^ 1, 
are logarithmically divergent at A; — 0, ^ 0, whereas the exact values 11(0, 0) 



and n(0, 0) are finite |jT3 



n(0,0) = -|^ = -^; 11(0,0) = -(4 -c^), (1) 



where n is the total concentration of bosons, /i the chemical potential, cb = ynVo/m 
the velocity of sound in the Bogolyubov approximation for a weakly nonideal Bose 
gas [0, Vq = V{0) the zeroth Fourier component of the potential, and c the speed 
of sound in the k — >■ limit for the spectrum of elementary excitations e{k) ~ |k|c 
in the Belyaev theory M: 



'Si2(0)/m* . (2) 

Here Si2(0) is the anomalous self-energy part of bosons at zero 4-momentum p = 
(k, e) = 0, and m* is the effective mass of quasiparticles, which is determined by 
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the relation 



1 

m* 



2 
B 



1 ^9Sn(0) 9Si2(0) 



2m 



d\k\ 



d\k\ 



(3) 



where Sii(O) and Si2(0) are, respectively, the normal and anomalous self-energy 
parts (at — > 0, e — > 0), and 



B 



1 



gSn(O) 
de 



-Sn(0) 



(4) 



The model of a dilute Bose system of hard spheres with a small parameter f3 = 

in which, by means of a summation of the 



n/kQ <^ 1, considered in Ref. | 
"ladder" diagrams with single-particle Green functions Go{p) in the zeroth approx- 
imation in f], it is possible to exclude the infinite repulsion, leads to a finite value 

of Si2(0), 

Si2(0) = , (5) 



m 



ao being the vacuum scattering amplitude of the particles and Uq the concentration 
of bosons in the BEG (po = mrio). 

At the same time, in Ref. [§, taking into account an exact thermodynamic 
equation 



5Sn(0) 
de 



drii 
. dm 



1 



no a/i 



(6) 



where rii = n — riQ is the concentration of supracondensate bosons, exact asymptotic 
relations 



G'ii(p^0) = -G'i2(p^0) 



nQiTic 



n{e^ — c^'k? + i5) 



(7) 

m an 



were obtained. 

However, it was shown in Refs. ||12|, ||13| (see also ||2^) that at p = (k, e) = the 
anomalous self-energy part is precisely equal to zero, Si2(0) = 0. Problems then 
emerge with the determination of the velocity of sound (Q) and the quasiparticle 
relaxation time r ~ S]^2^(0) [^, as well as with the asymptotic formulas for the 
normal Gii(p) and anomalous 6*12 (p) Green functions at p — > 



Gii(p^0) = -Gi2(p^0) 

because at Si2(0) = the relations (^ and I 

gSn(O) _ 
dt " ' 



Sl2(0) 



I reduce to the identities 



5 = 



(8) 



(9) 



so that Eqs. (^ and (|^) with account for (^ contain uncertainties of the 0/0 type. 

With the purpose of fixing these controversies, as well as the infrared divergences 
and nonanalyticities at p — emerging in the nonrenormalized theory, a renormal- 
ization procedure for the field perturbation theory was worked out in Ref. 
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employing the method of "combined variables" , which in the long- wave domain 
(|k| < ko) are in fact the hydrodynamic variables in the spirit of Landau quantum 
hydrodynamics while in the short-wave domain (|k| > ko) they reduce to the 



bosonic creation and annihilation operators. As was shown in Refs. |T^,[^, the 
perturbation theory built on such "adequate" field variables, does not suffer from 
infrared divergences at p — 0, whose source at T = is the divergence of long- wave 
quantum fluctuations (acoustic Goldstone oscillations), associated with a sponta- 
neous breakdown of continuous gauge and translational symmetries in the SF state 
of a Bose system with a uniform coherent condensate. Such oscillations are essen- 
tially the hydrodynamic first sound in liquid '^He, propagating with the velocity of 
ci ~ 236 m/s. 

The choice of combined variables 



28 



p9|] leads to the renormalized anomalous 
self-energy part Si2(p) which does not vanish at p = 0. Then one can formally 
restore all the results of the renormalized field theory 0,0, but this time in terms 
of the renormalized quantities Gik{p) and Sjfc(p), which do not contain singularities 
at p — > (save for the pole part Gikip) ~ bl~^)- In particular, the square of velocity 
of first sound Ci, in accordance with Ref. (0), at T — > must be equal to 



'dP\ _ Si2(0) 



m 



(10) 



where the derivative of the pressure P with respect to the total density p is taken 
at constant entropy a, and the renormalized effective mass m* is determined by the 
relations (^ and (|p with Sjfc(O) substituted for Sjfc(O). 

In view of the aforesaid, we will work with the combined variables 



where 



1 H ^t:^^ \- «0L 



2(nL) 



k|<fco 



ikr 



(12) 



{nL)e 



Such an approach means that the sepaartion of the Bose system into a macro- 
scopic coherent condensate and a gas of supracondensate excitations is made not 
on the statistical level, like in the case of a weakly nonideal Bose gas 0, but 
on the level of ab initio field operators, which are used to construct a microscopic 
theory of the Bose liquid. Note that the approximate expression for the long-wave 
part of the boson field operator \1/ is given with account only for first-order 
terms in the expansions over the slowly changing (hydrodynamic) phase 0l and the 
small deviation of the density nL from its mean value (^l). In Refs. [^], it was 
assumed that at low T, because of rather weak interaction (m/coV^(/co) 1), almost 
all the particles are in the Bose condensate, and therefore the value of momentum 
kg in Ref. |^ was chosen in such a way that the approximate equation (tzl) — uq 
(no being the concentration of particles in the BEG) take place. However, in a Bose 
liquid with strong interaction, when the single-particle BEG is strongly suppressed 
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[riQ <^ n), the value {hi,) should be normalized to the density Ug = Ps/m of the SF 
component. 

The system of Dyson-Belyaev equations which allows one to express the 

normal Gn and anomalous Gu renormalized single-particle boson Green functions 
in terms of the respective self-energy parts Sn and S12, has the form 



G'ii(p,e)= Go (-P, -e)-Si 



P, 



/Z{p,e); 



Here 



G'i2(p,e) = Si2(p,e)/Z(p,e) . 
-p, -e) - Sii(-p, -e) Go ^(P' e) - Sii(p, e) 



|Si2(p,e) 



G-' 



lP>ej 



2^ + ^" + '^ 



{S +0) 



(13) 
(14) 



(15) 
(16) 



where is the chemical potential of the quasiparticles, which satisfies the Hugengoltz- 
Peins relation |]^: 

/i = Sn(0,0)-Si2(0,0). (17) 

Due to a strong hybridization of the single-particle and collective branches of 
elementary excitations in the Bose liquid with a finite BEG (no 7^ 0), the poles of 
the two-particle and all many-particle Green functions, as well as the full four-pole 
function of the pair interaction of bosons r{pi,p2,P3,P4), coincide with the poles of 
the single-particle Green functions Gik{p,e) 0, 0. 

Therefore the spectrum of all elementary excitations with zero spirality is deter- 
mined by the zeros of the function Z[p, e): 



where 



^ + tl,{p,E{p))~p 



1 2 



1/2 



\f:,2ip,Eip))\H +s»i(p,E(p)), (18) 



^ [Sii(p,e) ± 



P, 



the (+) sign corresponding to the symmetric part of S^j^, the (— ) sign — to the 
antisymmetric part_S"^. In the sequel, we will assume that En is an even function 
of p and e, so that S^^^ = and S^^ = En. Then, relation (|l^) ensures the acoustic 
dispersion law for the quasiparticles at p ^ 0, and the first derivative of Sn over 
e in Eq. (D at e and p ^ vanishes: ^^iiM = q, so that B ^ 0, while 
^12(0) = TiQ-^, in accordance with Eq. (^. 

As was shown in Ref. |2^, for a Bose liquid with strong enough interaction 



between particles, when the BEG is strongly suppressed, one can, when defining 
Ejfc(p, e) in the form of a sequence of irreducible diagrams containing condensate 
lines [^], restrict oneself, with good precision, to the first (lowest) terms in the 
expansion over the small BEG density (no <C n). Such an approximation is exactly 
opposite to the Bogolyubov approximation for a weakly nonideal Bose gas with 
an intensive BEG, when no — n. 
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As a result, up to terms of first order in the small parameter riQ/n <^ 1, for a 
Bose liquid one gets the "trimmed" system of equations for Ej^ [EB]: 



Sn(p, e) = noA(p, e)V{p, e) + n^V{0) + ^n(p, e) ; 
Si2(p, e) = noA(p, e)V{p, e) + ^i2(p, e) , 



(19) 
(20) 



where 



—Gij{k)V{p - k, e - ^)r(p, e, k, a;) 



(21) 
(22) 



(27r)3 

r(p,e) = \/(p)[l-y(p)n(p,e)]-^ . 

Here V^(p) is the Fourier component of the input potential of pair interaction 
of bosons, V^(p, e) is the renormalized ("screened" due to multiparticle collective 
effects) Fourier component of the retarded (nonlocal) interaction; n(p, e) is the 
boson polarization operator: 

"^P'') = ^/(I^/^^(P'''^'") 
X {G'ii(k, cj)G'ii(k + p, e + cu) + G'i2(k, cu)G'i2(k + p, e + w)} 



(23) 



r(p, e; k, oj) is the vertex part (three-pole function), which describes multiparti- 
cle correlations; A(p,e) = r(p,e, 0,0) = r(0,0,p,e) , and rii is the number of 
supracondensate particles (ni ^ wq), which is determined from the condition of 
conservation of the total number of particles: 



n 



nQ + ni = riQ + i 



d^k r duj 



(2vr) 



27r 



Guik.u) 



(24) 



In the sequel, as well as in Ref. p6[, in the integral relations (|2T|) and (psf ) we will 
only take into account the residues at the poles of single-particle Green functions 
Gjj(p, e), neglecting the contributions of eventual poles of the functions r(p, e, k, u) 
and V^(p, e), which do not coincide with the poles of Gjj(p,e). As a result, taking 
into account relations (|l^)-(|l^), ([T8|) , ([T9| ) and (|^), Eqs. (|2T|) on the mass shell 
e = -E(p) assume the following form (at T = 0): 



^n(p,i5^(p)) = ^ 



(2vr) 



X V(p-k,i5;(p)-E(k)) 



3r(p,E(p);k,E(k)) 



A(k) 



E(k) 



(25) 



^12(P) =^12(P,^(P)) = -- 



d^k 



r(p, E(p); k, i?(k))\/(p - k, E(p) - E(k)) 



where 



(27r)3 

noA(k,E(k))y(k,E(k)) + ^i2(k) 
E(k) 



A(p) = noA(p, E{jp))V{p, i?(p)) + niK(O) + vl/n(p) + 1^ 



(26) 
(27) 
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Then the nonhnear equation ( [T8| ) for the quasiparticle spectrum E{p), according to 
Eqs. (ra), (PO), takes the form 



Eip) = JA^ip) - noA(p, Eip))Vip, E(p)) + ^^p] 



l2 



(28) 



and the total quasiparticle concentration in the Bose liquid is determined by the 
relation 



n 



(2^ 



(29) 



The Hugengoltz-Peins relation (|T^, according to Eqs. (|19]) and ([20|), can be 
represented as 

;i = niV(0) + ^n(0)-^i2(0) , (30) 
result of which Eq. (|2^ takes the form 

Aip) = noA(p, ^(p))V^(p, E{p)) + [§n(p) - ^11(0)] + ^12(0) + ^ . (31) 

From Eqs. (|28|) and (|3lD it follows that the quasiparticle spectrum, because of 
the analyticity of the functions \^jj(p, e), is acoustic at p — 0, and its structure at 
p 7^ depends essentially on the character of the renormalized pair interaction of 
bosons. 

If one assumes that the functions and \l/i2 depend weakly on the explicit 
form of the quasiparticle spectrum E{p), then by virtue of Eq. (^, one can, with 
good precision, assert B = 1, so that the effective mass, according to Eqs. (|]), ([T9|), 
(EO), is determined by the expression 



1 



1 

m 



92^11(0) 92^12(0) 



d\p\ 



d\p\ 



and the expression for the velocity of sound, according to Eqs. 
be cast in the form 

^12(0) 



A(0,0)\/(0,0)n/m* ; n = rio + 



(32) 

and (PH), can 
(33) 



A(0,0)V(0,0) 

which is analogous to the expression for the Bogolyubov velocity of sound for a 
weakly nonideal Bose gas cb '- 



'V{0)n/m. Then the conditions n > and c = Ci 
imply severe constraints on the choice of the parameters of the model of interaction 
of the bosons (see below). 



Indeed, at p ^ and E{p) = ci\p\ 0, Eq. (p6D, due to the the momentum 
dependence of the spectrum E{p) and the functions V"(p) = V^(p, -E(p)), A(p) = 
A{p, E{p)) and \E'i2(p) being isotropic, assumes the form 



^12(0) 



k'^dk 



(27r)2 J E{k) 



noA\k)V\k) + A(k)t^(k)^i2(k) 



(34) 



It follows that the first integral addend on the right-hand side of Eq. (^4|) is always 
negative, so that the value \^i2(0), which plays the role of the pair order parameter 
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can also be negative at not too small values of uq, regardless of the sign of 
the renormalized interaction A(k)V^(k). The condition \^i2(0) < means that the 
phase of the PCC is opposite to that of the BEC, because Uq > 0. Moreover, 
in this case, due to the condition A{0)V{0) > 0, which ensures that the system 
is globally stable against a spontaneous collapse, at sufficiently small densities of 
the BEC, in accordance with Ref. (|5BD, the values n and can become negative, 
which corresponds to an instability in the phonon spectrum, which occurs when 
1^12(0)1 >noA(0)V(0). 

However, if the pair interaction between bosons in a broad enough region of 
the momentum space has, due to some reason, the character of attraction, i.e., 
A(k)V^(k) < at k 7^ 0, and if the magnitude of that attraction is large enough 
(see below) and the BEC density is small enough (no -C n), the second (positive) 
addend on the right-hand side of Eq. (|3^ can outweigh the first (negative) one. 
Then \l/i2(0) will be positive, and the phase of the PCC will coincide with that of 
the BEC, so that h > and > 0. 

Note that when the BEC is totally absent (no = 0), when the integral equation 
(p6|) is homogeneous, i.e., degenerate with respect to the phase \^i2(p), the condition 
of stability of the phonon spectrum = \E'i2(0)/m* > is secured by the choice of 
the respective sign (phase) of the pair order parameter ^12(0) > (see [Q). 

On the other hand, since at T = the density of the SF component ps coincides 
with the total density p = mn of the Bose liquid, assuming n = n and with account 
for (|20|), one gets the following relations: 

Si2(0) 

Ps = PO + Ps = TTT- = ; jS) 

A(0)l^(0) ^ ^ 

^12(0) 

Ps = mni = m , [6b) 

^ A{0)V{0) ^ ' 

where po = mriQ is the single-particle BEC density, and ps is the density of the 
"Cooper" PCC. The concentration rii = n — uq is then determined from the rela- 
tion (P^D, and for hquid ^He at T — » 0, in accordance with the experimental data 
|jTT| , it should be not less than 90% of the full concentration of "^He atoms. Thus, 
the SF component in this model is a superposition of a single-particle and a pair 
coherent condensates, and relations (|29|) and (^) impose additional constraints on 
the parameters of the microscopic theory of the SF Bose liquid. 

3 Influence of a renormalized pair interaction on the spec- 
trum of elementary excitations in a Bose liquid with a 
suppressed BEC 

To describe interaction of helium atoms in real space, various phenomenological and 



semi-empirical potentials [^-[^ are conventionally used, which describe strong re- 
pulsion at small distances and weak van der Waals attraction at large distances (see 
Table I). As was shown in |J0| by means of a numerical solution of the Schrodinger 
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equation, some of those potentials give rise to bound states — discrete levels with 
very small binding energy {AE < 0.1 K). 

However, all these potentials are characterized by a divergence (±00) at r — 
and are not suitable for the description of pair interaction in momentum space, since 
the respective Fourier components are infinite. Moreover, in a condensed state one 
cannot directly use the free atomic potentials — it is rather necessary to address the 
problem of constructing adequate pseudopotentials [^. Therefore, as the input 



potential of pair interaction, a regularized "hard-spheres" potential |^3[, with 
the radius equal to the quantum radius of the helium atom, tq = 1.22 A, will be 
employed. 



As was shown in Refs. |2J], for spherically symmetric scattering (S'-wave), 
taking into account that the radial wave function u{r) of relative motion of particles 
vanishes at the infinite jump of the potential V{r) ^ 00 at r = 2.44 A, as well as 
the boundary condition V{r)u{r) = X6{r — a), the Fourier component of the pair 
potential is a finite sign-changing function of momentum transfer: 

Vip) = VoMpa) ; Ux) = ^ . (37) 

X 

Here jo{x) is the spherical Bessel function of zeroth order, and Vq is a positive con- 
stant which is determined in a self-consistent way from a nonlinear integral equation 
for the single-particle Green function at p — > and depends on the dimensionless 



Bose liquid density na (see Refs. [^, p^). 



The potential (|37|) is shown with a dashed line in Fig. 1. It corresponds to 
repulsion, V{p) > 0, in those regions of the momentum space where sin(pa) > 
(in particular, when pa < vr), and to attraction, V{p) < 0, in those ones where 
sin(pa) < (for example, n < pa < 2tt). 

As a different example, consider the Fourier component of a model potential in 
the form of a "Fermi" function 

nr) = K,|exp(^^:^^+l| , (38) 

which at 6 = degenerates into a "step" of a finite height Vq at r < a, corresponding 
to a model of semitransparent spheres [|25l. In this case the Fourier component is 
expressed in terms of the first order spherical Bessel function: 

Vip) = V,'^; ^^^^^ = Mx)-xcos{x) _ ^3g^ 
pa x'^ 

The Fourier component of a smooth potential V{r) in the form of a Lindhardt 
function [^], whose derivative at the inflection point r = a is —00, has the same 
shape (see Appendix A). The sign-changing oscillations of the Fourier component 
V{p) in the momentum space are then formally analogous to the Ruderman-Kittel- 
Ksui-Yoshida oscillations in exchange interaction of spins or to the Friedel 



oscillations of the screened Coulomb potential with period Tr/kp in real space, 
which arise as a result of scattering of electrons (fermions) on a Fermi sphere, filled 
according to the Pauli principle, of diameter 2kF {kp being the electron Fermi 
momentum) . 
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It should be stressed that the sign-changing Friedel oscillations are not directly 
associated with the jump of the Fermi electron distribution function at T = 0, 
but rather with a relatively weak singularity (inflection point) in the momentum 
dependence of the static polarization operator nc(p, 0), which is characterized by a 
logarithmic divergence at p = 2kp of the first derivative of nc(p, 0) with respect to 
p. This means that the sign-changing oscillations of the Fourier component of the 
potential V{r) can arise not only in the "hard spheres" model with an infinite jump 
of V{r), but also with finite jumps of the potential, or for "smooth" potentials with 
weak fractures or inflection points with an infinite derivative with respect to r. 

Thus, the existence of negative values V{p) < 0, i.e., of an effective attraction 
in some regions of the momentum space, is a consequence of certain peculiarities 
(jumps, fractures, inflection points) in the radial dependence of the pair interaction 
potentials, and is of quantum mechanical nature. 



If one substitutes the oscillating potential (|37|) (or (|39|)) into the Bogolyubov 
spectrum of a dilute quasi-ideal Bose gas pi 



Eb{p) 



p1 

2m 



|l + 2„y(p) 



1/2 

' , (40) 



then, by choosing two parameters, Vq and a, independently, one can achieve a rather 
satisfactory coincidence of the spectrum Eb{p) with the elementary excitation spec- 
trum Eexpip) in liquid ^He derived from neutron scattering experiments (Fig. 2, 
solid curve). However, the self-consistent solution for a = 2.2 A and na^ ~ 0.23, 
which was obtained in Refs. [^, |2^ in the framework of the regularized "hard 
spheres" model, differs considerably from ii^exp(p), and for a more realistic param- 
eter a = 2.44 A, characteristic of ^He, when na^ = 0.315, the spectrum (^) with 
the potential ( P7| ) turns out to be unstable, because -E'^(p) < in some range 
of p (Fig. 2, dashed curve), which points out that the Bogolyubov theory is 
unapplicable for the description of the "^He Bose liquid. 

On the other hand, multiparticle collective effects in the Bose liquid, according 
to Ref. (p2l), lead to an essential renormalization of the pair interaction, which 
determines the normal and anomalous self-energy parts, Eqs. ([I9D and (|20|) . Taking 
into account Eqs. (|2^) and (|37D , the vertex part (four-point function) of the retarded 
interaction between bosons assumes the form 

T>^ N Vosin(pa) 

V{p,uj) = — . . . (41) 

pa — VoLL[p, 00) sm(paj 

where n(p, u) is the bosonic polarization operator (^) in the SF state, which, with 



account for the pole parts of the Green functions ([T3|) and ([T^ , is calculated in 
Appendix B. 

An important feature of the renormalized interaction (^) is that in those regions 
of phase volume (p, cu) in which n(p, ti;) < 0, the repulsion (when sin(pa) > 0) gets 
suppressed while the attraction (when sin(pa) < 0) gets effectively enhanced. 

As follows from Eqs. (^) and (|3TD , the main influence on the quasiparticle 
spectrum -E'(p) comes from the shape of the interaction potential (^) on the "mass 
shell", when uj = E{p). The real part of the polarization operator then has the 
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form (see Appendix B): 

Ren(p,E(p)) 
f ^-(k,p) 



(27r)3 E(k) - E(k - p) - ^(p) 



(42) 



E{k) [E{k) + E{k-p)- E{p)] E{k) [E{k) + E{k-p) + E{p)] 



where 



i^-(k,p) 



E(k) + i^-/i + Sn(k,i5;(k)) 
2m 



X 



Eik) - Eip) + l^^l! _ ^ + Sn(k - p, Eik) - E{p)) 



2m 

Si2(k, E(k))Si2(k - p, E{k) - E{p)) 



(43) 



i^+(k,p) 



E{k-p) + 



2m 



-/i + Sn(k-p,E(k-p)) 



E(k - p) + E{p) + — - + Sn(k, E(k - p) + E{p)) 
2m 



(44) 



+ Si2(k - p, E(k - p))Ei2(k, E{k - p) + ^(p)) . 

The key feature of the screened potential V^(p, -E'(p)) is the fact that Re (n(p, E{j))) < 
for all p > 0, if the quasiparticle spectrum E{p) is stable with respect to decays 
into a pair of quasiparticles [0, i.e., if for all p and k the conditions 

E{p) <E{k)+E{k-p) ; E(k) <E{p) + E(k-p) (45) 

are fulfilled. 

Indeed, as follows from Eq. (0), the common denominator in front of the curly 
braces is always negative, 

[E{k)-E{k-p)-E{p)]<0, (46) 
whereas the denominator in the first term in the curly braces is always positive, 

[E{k)+E{k-p)-E{p)]>0 (47) 
and smaller than the positive denominator in the second term 

[E(k) + E(k - p) + E{p)] > . (48) 

At the same time, as numerical simulation has shown, the numerators F_ and 
in both terms remain positive for any p and k (see Fig. 3a). Therefore, the common 
sign of integrands in Eq. ( ^2]) is negative, so that lie {Il{p, E{p)) < 0. Figure 3b 
shows the p dependence of Re (n(p, E{p)), obtained numerically by iterations. 

Taking into account the negative sign and the relatively weak momentum depen- 
dence of Re (n(p, E{p)) in a wide region p ^ 0, one can approximate the renormal- 



ized potential (|41| ) at cu = E(p) in the integrands of Eqs. (|25D, (p6D with a more 
simple one: 

- . ^-'Pf (49) 

pa + a sm[pa) 
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where a = Vo|n|, and |n| = |Re (n(p, E{p))\ is the mean absolute value of the real 
part of the polarization operator on the mass shell in the domain of existence of 
the spectrum E{p). In what follows, all the numerical simulations are based on 
the model potential in which the quantities Vq and a are free fitting parame- 
ters. Figure 1 depicts the function V{p)/Vo for different values of the dimensionless 
parameter a. As far as the vertices A and F are concerned, their relatively weak 
p and u dependence may be disregarded, putting A ^ F ~ A(0, 0) = const and 
incorporating the constant value of A(0, 0) into Vq. 



As a result, 
form. 



equations (pSD and (p6D for the functions boil down to a simpler 



^11 (P) 



V[p 



(2vr) 



k) 



AM 
E(k) 



1 



^12(P) 



V {p — K) 



(27r) 



i?(k) 



where 



E{p) = \ Al{p)- noV^(p) + M/i2(p 



Ao(p) = noV{p) + 
In this case, for p — 0, from Eq. 



^ll(p) -^ll(0)J +^12(0) + 

taking into account 



2m 
and 



E{p 0) = IpK/ noV^(O) + ^12(0) /m* 



(50) 
(51) 

(52) 

(53) 
, one gets 

(54) 



where V^(0) = Vo/[l — Von(0,0)], and 11(0,0) is determined by the first relation in 
(ffl), so that for the velocity of sound one has the following relation (see (|55|)): 



1 

fh* 



UoVq 



1 + noVo/mc^ 



+ ^12(0) 



(55) 



which is equivalent to a biquadratic equation for c. Under the conditions nVo 3> 



mc 



no <^ n, from Eq. (55) there follows an approximate expression 



^12(0) 



ra 



(56) 



The parameters Vq and a were chosen in such a way that the phase velocity E{\p\ 
0)/|p| coincide with the hydrodynamic sound velocity ci ~ 236 m/s in liquid ^He. 
On the other hand, the choice of those parameters was to ensure the maximal 
coincidence of the spectrum E{p) with the experimental spectrum ii^exp(p) in ^He 
I^-|T^ at all values of p. 

Figure 4 (a,b,c) depicts the momentum dependence of the functions \^ii(p), 
^12 (p) and Ao{p), calculated according to Eqs. ([5(]|), (^Tp and ( p5D for the values of 
a = 4.4, Vo/(47r^a^) = 17.7 K at a = 2.44 A, and Fig. 5a presents the quasiparticle 
spectrum E{p) obtained from Eq. (^) for those same values of the parameters. It 
can be seen that this spectrum is in qualitative agreement with the experimental 
spectrum E^^^p in ^He, however the numerical correspondence of the positions and 
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values of the maximum and minimum of quasiparticle energy, -Emax = 15.9 K at 
Pmax = 1.05 A and -Emin = 10.7 K at pmin = 1.72 A , cannot be deemed sat- 
isfactory. Besides, the sound velocity calculated according to Eqs. (|33D and ( ^5l) 
turns out to be too low, c = 2.08 x lO'* cm/s, while the total particle concentration 
calculated according to the formula (^91) is too high, n = 2.57 x 10^^ cm~^, the 
concentration of the particles in the BEC being low, no = 0.03n. 

There is a much better numerical agreement at the values of a = 4.52, Vq/ (47r^a^) = 
11.2 K and a = 3.0 A. The respective spectrum E{p) is shown in Fig. 5b and 

is characterized by the following values: -Emax = 14.28 K at pmax = 1.13 A , 
Emin = 9.8 K at Pmin = 1.95 A c = 2.34 X 10'^ cm/s, the total concentration be- 
ing n = 2.18 X 10^^ cm~^ and the BEC concentration uq = 0.09n (which agrees with 
the experimental data [|ll|]). Such an agreement should be deemed quite satisfactory, 
taking into account the simplified model of pair interaction employed. 

It is seen from Fig. 4-5 that the nonmonotonous nature of the spectrum E{p), 
including, in particular, the presence of the "roton" minimum, is determined mainly 
by the momentum dependences of the functions \Efii(p) and Aq{p), which have deep 
minima due to the oscillations of the sign-changing potential V{p) in the domain 
p < 2iT/a (see Fig. 1). The theoretical spectra obtained are in good agreement with 
the experimental spectrum of ^He both in the positions and the absolute values of 
the maximum and minimum of E{p). 

One should stress that at comparatively small values of the parameter a it is 
impossible to render the quantity ^1^12(0) positive, because repulsion between bosons 
prevails over effective attraction when integrating over p in Eq. (0). At the same 
time, because of the BEC density being small {uq <^ n), the quantity Si2(0) = 
r2oA(0)V^(0) — |\t'i2(0)| becomes negative, which corresponds, according to (|TD|) and 
(P3|), to a phonon instability of the spectrum (c^ < 0). 

On the other hand, at large enough values of a, the effective attraction V{p) < 
in the region tt < p/a < 2tt (Fig. 1) turns out to be so strong that around the 
negative minimum of V{p), the radicand in Eqs. (pSf ) and (^) becomes negative, 
i.e., the quasiparticle spectrum E{p) becomes absolutely unstable (imaginary) in 
the region in question, analogously to the Bogolyubov spectrum with a nonrenor- 
malized potential in the "hard spheres" model (see Fig. 2). Therefore, the domain 
of parameters Vq and a in which there exists a quasiparticle spectrum (^) which 
would be stable at all p and be in agreement with experiment, is rather narrow. 
It is possible, in principle, to consider the inverse problem and find the effective 
"pseudopotential" of pair interaction of particles in the "^He Bose liquid, starting 
from the shape of the empirical quasiparticle spectrum Ee^pip), however this is quite 
tedious. 

4 The criterion of superfluidity and limiting critical veloc- 
ities 

We conclude with a brief discussion of the applicability of the Landau criterion of 
superfluidity to He II and of the value of limiting critical velocity in the absence of 
quantum vortices in the Bose liquid where the BEC and PCC coexist. 
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As was noted in Introduction, the elementary excitation spectrum £^exp(p) as 
observed in neutron scattering experiments |]T3|-|]TB|, leads to a rather overestimated 
value of the critical velocity, determined by the roton minimum (in accordance with 
the Landau criterion), as compared with the experimentally measured velocities of 
destruction of the SF flow. This has to do with the emergence of quantum vortices 
and vortex rings in He II IH], however under conditions when creation and/or 



motion of those vortices is hindered, the critical velocities increase sharply pU[ , 
21 1, and at low temperatures < IK they can assume values comparable with 
Vc = min [e{p)/p] |22| . 



One should stress that this situation is quite analogous to the one with the type 
I superconductors, in which the critical current jc is determined by the condition 
of creation and pinning of Abrikosov quantum vortices on the surface of the super- 
conductor or near various defects of the crystal lattice, whereas the true maximal 
value of jc — the so-called decoupling critical current, which is determined by the 
process of decay of Cooper pairs |^ , — is much larger and can only be observed in 
thin wires, whose thickness is much less than the London depth of penetration of 
magnetic field into the superconductor, which hinders vortex creation. 

In the ^He Bose liquid at finite temperatures, T 7^ 0, together with the spectrum 
Eexpip), which at p ^ corresponds to the first (hydrodynamic) sound with phase 
velocity ci, there arises — due to the emergence of the normal component p„ in 



He II — the second sound, whose velocity C2 <C ci at T > 1 K. It is known |1I0 



^ that, due to heat expansion of liquid ^He being weak, the second sound branch 
reduces basically to oscillations of temperature (entropy) without any noticeable 
transfer of total mass of the normal (p„) and SF (p^) components, whose oscillations 
are in counterphase. Therefore, the excitations of second sound with energy e2(p) = 
C2P cannot be observed in standard neutron scattering experiments, unlike the first 
sound ei{p) = Cip, which corresponds to cophase oscillations of the densities ps 
and Pn- At the same time, in the two-component Bose liquid there may coexist 
two different types of acoustic Goldstone excitations, associated, on the one hand, 
with spontaneous breakdown of gauge symmetry due to the degeneracy over the 
phase of the coherent SF condensate at T — > 0, and the breakdown of continuous 
translational symmetry, i.e., of the homogeneity of the total density p = Pn+Ps (first 
sound); on the other hand, with spatially inhomogeneous deviations of temperature 
T from the uniform distribution due to the oscillations of the density of the gas 
of normal excitations (second sound). Nevertheless, the second sound excitations 
do carry over some energy and therefore have to be taken into account in the 
determination of the minimum critical velocity, according to the original concept of 
the Landau supefiuidity criterion [|l|], which includes all the types of excitations in 
the quantum liquid. 

In this respect, one can assume that at those temperatures where C2(T) < Vc — 
60 m/s, the maximal allowed critical velocity of the macroscopic SF flow in He II in 
the absence (or with strong pinning) of quantum vortices cannot exceed a number of 
the order of the second sound velocity C2(T), which tends to zero as T — > Tx, together 
with the SF component density Ps(T). Precisely such a situation is characteristic 
of superconductors, in which the decoupling critical current turns into zero in the 
critical point T = Tc together with the energy gap A in the quasiparticle spectrum. 
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Finally, it is worth noting that the coexistence of a weak BEC and an intensive 
PCC conserves the integer value of the quantum of circulation of the SF velocity 
in the vortices n = h/m, due to the total mutual coherence of those condensates 
in the SF component ps- Indeed, a sufficiently strong effective attraction for the 
screened Fourier component of the singular "hard spheres" potential provides for 
the formation of a condensate of bound bosonic pairs with a positive sign of the 
pair order parameter \E'i2(0), whose phase in that case coincides with the one of the 
BEC. 



5 Conclusions 

Thus, employing the renormalized field theory for the description of the SF state 
of the Bose liquid at T — > with account for a small density of the single-particle 
BEC allows one to formulate a self-consistent model of superfiuidity, in which the 
SF component at T ^ is a coherent superposition of the single-particle BEC, 
suppressed due to interaction, and an intensive PCC, which arises due to an effective 
attraction between the bosons in the momentum space. Such an approach lets one 
obtain, in the framework of the "hard spheres" model, an explicit form of the 
quasiparticle spectrum, which, with a suitable choice of parameters, coincides with 
good precision with the experimental spectrum of elementary excitations in ^He. 
We are sincerely grateful to P.I. Fomin for numerous useful discussions. 



6 Appendix A 



In order to calculate the Fourier components of the empirical potentials presented 
in the Table, it is necessary, in order to avoid the divergence at r — 0, to make a 
cutoff at some distance Tc, representing the diameter of a "hard core" . For example, 
for the Lennard- Jones potential, up to the main terms, one has 



oo 

U{p) = Ane J rdr sm(pr) 
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where Si(a;) is the integral sine. This potential can be extended to the region r < 
of the regularized Fourier component of the infinite repulsion potential in the "hard 
spheres" model |^ , or of the finite potential of "semi-transparent" spheres |^ . 
Such an approach is analogous to the method of pseudopotential in solid state theory 
W^ - Sign-changing oscillations of the Fourier component (A.l) are associated with 
a sharp fracture and a finite jump of the potential at the point r = r^. In particular. 
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Eq. (A.l) by itself corresponds, in coordinate space, to the model of pseudopotential 
with an "empty skeleton", when the potential is U{r) = aX r < r^, and adding 
the Fourier components of the potentials of "hard" or "semitransparent" spheres to 
that equation corresponds to the model of "filled skeleton" . 

Oscillating sign-changing Fourier components are also characteristic of certain 
types of smooth potentials with inflection points. As an example, apart from the 
Fermi function (PBI), there is a potential in the form of a static Lindhardt function 



2% in real space: 



W{r) 



Wo 
2 



l+x 



X 



X = r/ro 



(A2) 



which has an inflection point with an infinite negative derivative at r = tq. 
Fourier component of this function can be calculated exactly, and equals 



The 



W{p) = 



pro 



(A3) 



i.e., its momentum dependence coincides with the one of the Fourier component of 
the potential of semi-transparent spheres with a finite jump at r = Tq (see (P^). 
This is formally equivalent to the calculation of the Ruderman-Kittel-Ksui-Yoshida 
static oscillations |^| in real space for indirect exchange interaction of spin in metals 
with the Fermi distribution function at T = and with the Fourier component of 
spin susceptibility of electrons in the form of Eq. (A. 2) with x = p/ {^kp) [kp being 
the Fourier momentum of the electrons). 

At the same time, Fourier components of smooth continuous potentials with finite 
derivatives, nonsingular at r = 0, are characterized by weak oscillations. Consider 
for example the modified nonsingular Buckingham potential (see Table), in which 
the divergent term (— 1.49/r^) is replaced, at r < tq = 2.61 A, with an exponential 
{—Be~^^), finite at r = 0. The parameters, B = 0.57 eV and (3 = 1.84 A 
chosen in such a way that at r < tq, the expression 



are 



VUr) = Ae 
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i -1 



with A = 481 eV, a = 4.6 A , would vanish at r = tq, and its derivative over r at 
that point would be equal to that of the potential in the r > tq range: 
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at C = 611 eV, D 
component of this potential is 
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where 
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Evidently, the oscillation amplitudes in this case are exponentially small 
pared to the smooth part of the potential. 



(A8) 
as com- 



7 Appendix B 

The polarization operator (p3D can be calculated without account for the vertex 
part r, making use of the expressions ([T3|)-(p!6D in the form 



/d k 
^[/n(p,k,u;) + /i2(p,k,o;)] , (5.1) 

where 



/dz ~ 
— Gij(k, z)Gij{^^ - p, 2; - cj) . (5.2) 
ZTT 

Assume that the Green functions Gij have only one pole within the integration 
contour and are equal to 

Calculating the integrals (B2) with account for the poles at the points e = E{k) 
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and e = E{k — p) + uj yields 
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In the statistical limit (cj — > 0, p 0), expression (B5) reduces to 
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It follows that in a large region of momentum space, Jii(0,k, 0) < 0. The same 
result is obtained for the function (B.6) at p = and u = 0, i.e., /i2(0, k, 0) < 0, so 
that the static bosonic polarization operator 11(0, 0) is negative, which corresponds 
to a suppression of the "screened" repulsion at p — » 0. From Eqs. (B.5) and (B.6) it 
can also be seen that on the mass shell u = E{p), the integrals In and /12 remain 
negative in a wide region of momentum space because of the negative sign of the 
common denominator [E{k) — E{k — p) — E{p)] < and the positive sign of the 
denominator [-E'(k) + £^(k — p) — E{p)] > due to the fact that the quasiparticle 
spectrum E{p) is decayless [conditions (^5])]. 

Thus, the polarization operator on the mass shell e = E{p) is negative on the 
whole range p < 27r/a, where sin(pa) < 0. It is necessary to stress that this feature 
of the polarization operator, n(p, E{p)) < 0, is only characteristic of Bose systems, 
in which the single-particle and collective spectra coincide with each other and are 
counted from the common zero of energy, unlike the Fermi systems, in which the 
single-particle excitation spectrum begins at the Fermi energy, due to the Pauli 
principle. Therefore, in the Fermi liquid (^He) there can be no corresponding effec- 
tive enhancement of the negative values of the "input" interaction potential V{p), 
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so that the formation of Cooper pairs is only possible for nonzero orbital momenta, 
due to the true weak van der Waals attraction between fermions. Apparently, it is 
this fact that has to do with the critical temperatures of the SF transition in ^He i 
^He differing by three orders of magnitude. 
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Figure Captions 

V(pa) 

Fig. 1. The dependence of the ratio — — — on the dimensionless momentum 

ii for different values of the dimensionless parameter a = Vo|n|: Dashed curve 1, 
the "hard spheres" input potential at a=0; sohd curve 2 corresponds to a=2, curve 
3 — to a=3] curve 4 — to a=3.5. 

Fig. 2. The Bogolyubov spectrum (^) for a dilute quasi-ideal Bose gas, obtained 
by substituting the potential (^) with an independent choice of the two parameters 
a = 2.5 and Vo/(47ra^) = 169.9 K at a = 2.44 A. The solid curve corresponds to 
the instability of the Bogolyubov spectrum obtained by substituting the potential 
(0) into (|0|) for = 2.17 x 10^^ cm'^ i a = 2.44 A, characteristic of ^He. 

Fig. 3. (a) The numerators -F_|_ i F_ of both terms in expressions (^) and ( ^41) 
remain positive at any p and k. (b) The p dependence of the polarization operator 
n(p,£'(p)) ( ^21) on the mass shell, obtained numerically by iterations. 

Fig. 4. The momentum dependence of the functions (a) ^ii(p), (b) "^u^p), 
(c) Ao{p), obtained according to Eqs. (^), (0) i (jS^), for Vo/{47T^a^) = 17.7 K, 
a = 4.4. 

Fig. 5. (a) The elementary excitation spectrum E{p) obtained from Eq. (|5BD, for 
the same set of parameters as on Fig. 4. The maximum energy in the quasiparticle 
spectrum (|5^ a,t p/fi = 1.05 A is -Emax = 15.9 K, and the minimum one, at 
the roton minimum at p/h = 1.72 A , is i^min = 10.7 K. The hydrodynamic 
sound velocity Ci = 2.08 x 10^ cm/s and the total quasiparticle concentration n = 
2.57x 10^^1/cm^ at the BEC concentration Hq = 0.03n and the diameter a = 2.44 A, 
found according to Eqs. ( |55| ) and (b) The spectrum E{p) for Vo/(47r^a^) = 

11.2 K, a = 4.52. The maximum energy of the spectrum ( p2| ) at p/h = 1.23 A 
is -Emax = 15.85 K, and the minimum one, at the roton minimum ai p/fi = 2 A 
is Emin = 9.81 K. The hydrodynamic sound velocity Ci = 2.34 x 10^ cm/s and the 
total quasiparticle concentration n = 2.18 x lO^^l/ cm^ at the BEC concentration 
no = 0.09n and the diameter a = 2.95 A, found according to(|33|) and p^). The 
dashed line corresponds to the experimental curve. 
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